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Abstract 


Rotordynamic analyses have typically considered 
only the rotating assembly and in the last 10 years 
have been expanded to include housing dynamics. Ad™ 
vanced rotating machinery designs conceptually define 
the rotor well ahead of the housing. Preliminary 
rotordynamic evaluations are performed with a housing 
stiffness assumption that is typically determined only 
after the hardware is built. In addressing rotor sta- 
bility, a rigid housing assumption has shown to predict 
an instability at a lower spin speed than a comparable 
flexible housing analysis. This rigid housing assump- 
tion therefore provides a conservative estimate of the 
stability threshold speed. A flexible housing appears 
to act as an energy absorber and dissipates some of 
the destabilizing force. The fact that a flexible 
housing is usually asymmetric and considerably heavier 
than the rotor has been related to this apparent in- 
crease in rotor stability. Rigid housing analysis is 
proposed as a valuable screening criteria and may save 
a significant amount of time and money in construction 
of elaborate housing finite element models for linear 
stability analyses. 


Introduction 


A primary consideration for high-pressure turbo- 
machinery such as used in liquid rocket engines is 
minimum weight. To achieve this while delivering the 
pressures and flows required by the system, the design 
trend is toward small diameters and very high speeds. 
The results of this trend are rather spectacular as 
shown in Fig. 1. As an example, the Space Shuttle 
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Fig. 1. Turbopump Power Growth 


Main Engine High-Pressure Fuel Turbopump (SSME HPFTP) 
is one of the most advanced machines in the world with 
a power density of 100 horsepower/pound. , While there 
is a substantial weight savings in the rotor with its 
small diameter , much of the savings is associated with 
the lightweight housing. A lightweight housing is de- 
signed primarily as a pressure shell with a minimum of 
additional structure added to support the rotor. At 
high speeds, housing motion becomes a significant fea- 
ture affecting rotor critical speeds, stability, and 
bearing loads. 

Housing designs are dictated by pressure, flow, 
hydrodynamic efficiency, structural loads, and assem- 
bly requirements. While low speed pumps with massive 
and stiff housings are generally used to discuss rota- 
ting machinery problems, advanced lightweight housings 
are not as simply treated. The flexibility of the 
housing, coupled through bearings and seals, results 
in a considerably more complex problem, which is not 
as amenable to simple rules of design as in the past. 
Rotor stability is much more difficult to accurately 
predict today due to the added variable of housing 
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dynamics and the interaction with the coupling mechan- 
isms. These elements affect critical speeds and thus 
the rotor stability for operation above the rotor cri- 
tical speed. The concern is with the effect of housing 
flexibility on the stability analysis, or can more con- 
servative results be obtained in some other fashion for 
advanced rotating machinery during the preliminary 
study phase. The intent of this paper is to investi- 
gate the influence of housing flexibility on the rotor 
stability. 


Approach 

Initially, the housing dynamics influence on sev- 
eral of Rocketdyne rotordynamic analyses was determined . 
To accomplish this, compatible stability cases were re- 
peated with a rigid or ground housing assumption to 
quantify the change in threshold speed with elimination 
of housing flexibility and dynamics. A total of six 
different turbopump designs were included in this sur- 
vey including four SSME designs and two advanced rocket 
engine designs . The specific turbopumps included were: 

1. High-pressure oxygen turbopump (HPOTP) of the 
SSME 

2. Uprated high-pressure oxygen turbopump of the 
SSME for full power operation, referred to as 
the HPOTP Phase II 

3. High-pressure fuel (liquid hydrogen) turbopump 
(HPFTP) of the SSME 

4. MK38-0 redesign configuration that never mate- 
rialized for a modified SSME 

5. MK48-F high pressure fuel (liquid hydrogen) 
turbopump for the Orbital Transfer Vehicle 
(OTV) 

6. MK49-F high pressure fuel (liquid hydrogen) 
turbopump for the Advanced Space Engine (ASE) 

These six turbopump analyses were chosen because each 
has a detailed finite element model of its housing . The 
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Fig. 2. High-Pressure Oxygen Turbopump 


following is a brief description of each of the turbo- 
pumps included in the survey. 


The HPOTP is shown in Fig. 2 and a summary of its 
performance characteristics is provided in Table 1. 
Basically, the HPOTP consists of a double-entry cen- 
trifugal impeller supplying liquid oxygen to the main 
chamber. A portion of this discharge flow is diverted 
to a smaller centrifugal impeller mounted on the same 
shaft that supplies liquid oxygen to the preburners. 

A two-stage impulse turbine, which overhangs the bear- 
ings, drives the pump to a maximum operating speed of 
27,900 rpm. Two pairs of duplex angular contact ball 
bearings support the rotor. Preload springs between 
the bearing outer races in each pair maintain the axial 
load independent of rotor axial position. The turbine 
hot gases are isolated from the pump by a series of 
floating ring seals. Axial thrust on the rotor is 
reacted by a double-balance piston built into the main 
impeller. 
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Table 1. SSME HPOTP PERFORMANCE DATA 


KEY PERFORMANCE PARAMETERS BASELINE DATA JAN. 1984 


RPL 

FPL . 

MAIN 

BOOST 

MAIN 

BOOST 

PUMP INLET FLOWRATE, LB /SEC 

1070.6 

111,6 

1157.8 

129.4 

PUMP INLET PRESS, PSIA 

379.3 

3985.9 

392,4 

4403.6 

PUMP DISCHARGE PRESSURE, PSIA 

4108.7 

7106.7 

4556.0 

7861.2 

PUMP EFFICIENCY 

0.684 

0.809 

0,650 

0.800 

TURBINE FLOWRATE, LB/SEC 

61 

.8 

69.0 

TURBINE INLET PRESSURE, PSIA 

5015.3 

5660.8 

TURBINE INLET TEMP, R 

1407.2 

1596,3 

TURBINE PRESS RATIO 

1.506 

1.550 | 

TURBINE EFFICIENCY 

0.749 

0.755 | 

TURBINE SPEED, RPM 

27102 

29675 

TURBINE HORSEPOWER 

22902 

29174 


The uprated HPOTP Phase II design is basically the 
same as the original HPOTP of Fig. 2 except for a stif- 
fened shaft and straight annular seals with a smooth 
rotor and roughened stator at the small impeller wear 
rings . These modifications permit full power level and 
extend the operating speed to 30,000 rpm. 

The HPFTP is shown in Fig. 3 and a summary of its 
performance characteristics is provided in Table 2. 

The HPFTP consists of a three-stage centrifugal (liquid 
hydrogen) pump driven by a two-stage turbine to deliver 
the engine coolant. The turbopump operates at a maxi- 
mum speed of approximately 37,000 rpm. Two pairs of 
duplex angular contact ball bearings support the rotor 
radially and a thrust bearing is provided for start and 
shutdown axial positioning. Rotor masses are located 
inboard of the bearings . A balance piston built into 
the third stage impeller reacts the axial thrust during 
steady-state operation. Straight smooth annular seals 
between the impeller stages ensure the rotor stability. 

The MK38-0 redesign configuration included was a 
proposed concept for upgrading of the SSME to increase 
its thrust rating. Basically, this design has a direct 
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Fig . 3. SSME High-Pressure Fuel Turbopump 
Table 2. SSME HPFTP Performance Data 


KEY PERFORMANCE PARAMETERS 
1 PUMP INLET FLOWRATE, LB /SEC 

RPL 

146.4 

FPL 

1617 

PUMP INLET PRESSURE, PSIA 

204.1 

214.0 

PUMP DISCHARGE PRESSURE, PSIA 

6254.8 

7038.8 

PUMP EFFICIENCY 

0.758 

0.757 

TURBINE FLOWRATE, LB /SEC 

147.5 

164.0 

TURBINE INLET TEMPERATURE, R 

RvSlll 

1989.2 

TURBINE PRESSURE RATIO 

1.522 

1.558 

TURBINE EFFICIENCY 

0.770 

0.780 

TURBINE SPEED, RPM 

34,931 

37,076 

1 TURBINE HORSEPOWER 

63,288 

77,142 


drive, two-stage turbine and two centrifugal impellers 
mounted on a shaft supported by angular contact ball 
bearings . The maximum operating speed is approximately 
26,000 rpm. 
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The MK48-F turbopump is shown in Fig. 4 and a 



summary of its performance characteristics is provided 
in Table 3. This turbopump consists of a three-stage 
centrifugal (liquid hydrogen) pump driven by a two- 
stage overhung turbine. The pump design is similar to 
the HPFTP on a much smaller scale and has a maximum 
operating speed of approximately 95,000 rpm. Again, 
two pairs of angular contact ball bearings support the 
rotor radially and a balance piston built into the 
third stage impeller reacts the axial thrust. 


4. MK48-F Liquid Hydrogen Turbopump 
Table 3. MK48-F Performance Data 


KEY PERFORMANCE CHARACTERISTICS 


FLOWRATE, GPM 


DISCHARGE PRESSURE, PSI 


DESIGN SPEED, RPM 


TURBINE POWER, HP 


WEIGHT, POUNDS 






The MK49-F turbopump is shown in Fig. 5 with a 
performance characteristic summary given. This turbo- 
pump is very similar to the MK48-F but on a slightly 
smaller scale. The MK49-F has a maximum operating 
speed of 110,000 rpm and relies on straight smooth 
interstage seals to ensure the rotor stability. 



Fig. 5. MK49-F Liquid Hydrogen Turbopump 


Mathematical Models Description 

Existing finite element models of the rotor and 
housing for each of the six turbopumps included in the 
survey were used. The rotors are constructed of beam 
finite elements and the housings are typically composed 
of plate elements to include shell motion. Typical rotor 
and housing finite element models are shown in Fig. 6 and 
7, respectively, for the SSME HPOTP. The coupled sys- 
tem was analyzed using modal synthesis techniques that 
combine substructures of the rotor and housing-to- 
engine support structure connected through linear bear- 
ings, seals, and hydrodynamic interactions. This ap- 
proach has been described by Childs (Ref. 1) for rotor- 
dynamics applications. Generally, the housing-to-engine 
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NODAL PLACEMENT 

Fig. 6. SSME HPOTP Rotating Assembly 
Finite Element Model 



Fig. 7. SSME HPOTP Housing Finite 
Element Model 

support structure is unsymmetric and is lightly damped. 
These linear models not only include bearing and seal 
coupling but also, impeller diffuser interaction, tur- 
bine aerodynamic cross coupling, bearing damping, and 
rotor and housing structural damping. Output capabil- 
ity includes damped critical speeds and stability. 

Stability Calculations 

The Rocketdyne internal stability criterion is 
similar to that of Lund (Ref. 2) and Bansal (Ref. 3), 
which considers a mechanical system at steady-state 
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which undergoes a perturbation. The motion takes one 
of three forms. When the deviation from the original 
motion decreases with time, the system is stable. If 
the deviation increases with time, the system is said 
to be unstable. Neutral stability occurs when the 
motion of the system after a perturbation is oscilla- 
tory with no change in amplitude with time. 

The equation of motion for a typical mechanical 
system takes the following form: 

m y = -ky - cy + F(t) (1) 

Assume the solution to the homogeneous differential 
equation is as follows: 


~ Xt 

y = Qe 

substituting into the equation of motion: 


( 2 ) 


mX + cX + k = 0 

where the eigenvalue is: 


X = 


± Vc- 


4km 


2m 

X = -a ± bi 


where: 


2m 


b = 


\ Xkm - 


2m 


(3) 

(4) 

(5) 


Writing the general solution to the homogeneous equa- 
tion of motion: 


y = e ( 3 + bi)t + Q 2 e (_a “ bl)t (6) 


y = e 


K 


’ at e Mt + Q 2 e‘ blt ) 


(7) 
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Using Euler's formula, equation (7) can be rewritten 
in the general solution form: 

y = e at [Acos(bt) + Bsin(bt)] (8) 

Plotting the solution: 



Fig. 8. Rate of Amplitude Change (a > 0) 



t = TIME 


Fig. 9. Rate of Amplitude Change (a < 0) 


Figure 8 displays decreasing amplitude with time, and 
therefore a stable system for a > 0. However, when 
the real part of X is greater than zero, the system is 
unstable. This unstable motion is represented in Fig. 
9. Lastly, if a = 0, the e -at term becomes unity and 
the amplitude remains constant with time, as shown in 
Fig. 10. Therefore, the sign of the real part of the 
eigenvalue determines the stability of the system. 


a = 0 


TIME AFTER 
PERTURBATION 


Fig. 10. Rate of Amplitude Change (a = 0) 
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Stability Results of Survey 


As previously mentioned, tie stability results of 
a rigid versus a flexible housing assumption were com- 
pared, Figures 11 through 17 present rotor stability 



Fig. 11. HPOTP Stability Map of Turbine Overhang 
Mode (Hot Gas Seals Free to Float) 



Fig. 12. HPOTP (104%) Rotor Stability Map 
(Second Rotor Mode) 
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SPEED 


Fig. 13. HP OTP (109%) Rotor Stability Map 
(Second Rotor Mode) 



Fig. 14. MK38-0 (Redesign) Rotor Stability Map 

(1st Rotor Mode) 


maps for each of the six turbopumps described for only 
the potentially unstable inodes. The real part of the 
complex eigenvalue is plotted as a function of the 
rotor spin speed where a positive real part eigenvalue 
indicates the potential for unstable rotor motion to 
occur. 
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OPERATING 

SPEED 


Fig. 15. H: 



Figure 11 compares the lowest rotor mode stability 
of the HPOTP , which is the turbine overhang mode. With 
the hot gas seals free to float, the operating speed is 
far removed from the stability threshold speed for both 
a rigid and a flexible housing. A comparison of the 
second rotor mode stability of the HPOTP is shown in 
Fig. 12 and 13, which is the main impeller mode . Fig- 
ure 14 compares the turbine overhang mode for the MK38-0 









Fig. 17. MK48-F Rotor Stability Map (1st Rotor 
Mode - Floating Ring Seal Inactive) 


redesign configuration. A comparison of the lowest 
rotor mode stability was performed for the HPFTP as 
shown in Fig. 15. Figure 16 compares the first rotor 
mode stability of the MK49-F. With the turbine 
floating ring seals present, this turbine overhang 
rotor mode remains stable. The floating ring seal 
was deactivated to produce an unstable rotor mode for 
comparison purposes. A similar situation was found 
for the MK48-F. Figure 17 presents the. stability 
comparison of the first rotor mode, which is also a 
turbine overhang mode, with the floating ring seal 
Inactive. 

The stability maps described above indicate that 
as the rotor spin speed increases, the rotor in a 
flexible housing tends to approach an unstable condi- 
tion at a higher speed than the rotor in a rigid hous- 
ing. This observation is true for all six turbopumps. 
Thus, results of the survey to determine the effect of 
housing flexibility on rotor stability indicate that a 
rigid housing assumption is more conservative since 
it consistently predicts lower threshold speeds. It 
appears as though the housing acts as an energy ab- 
sorber. This can be seen when it is considered that 
the housing adds relative motion terms. 
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Consider the following: 


q = Acosoot + Bsinoit 
R 

q = Ccoscot + Dsinoot 

The elastic coupling force between the rotor and sup- 
port is: 

F = - (kq + cq) = -[k(q R - q g ) + C(q R - q g )] 

The work, or energy, generated per cycle is: 

" ‘C* ^ [k<,R " qs> + c<q R ' V 1 • 

Rearranging the above equation as follows: 

>-/ i- ( m r + c V] • % 

u 

2tt 

[-(kq s + cq s )].q R ^^-- 

Therefore it follows that if the housing support was 
rigid or ground, then q^ and q are zero. All the 
constant available energy in a S conservative system is 
then transmitted to the rotor and none dissipated in 
the housing. It follows then that a flexible housing 
absorbs some of the energy that would normally be 
present in the rotor. 

Another aspect of housing construction that aids 
stability is its typical asymmetry. Several technical 
papers discuss bearing support asymmetry related to 
stability gain. Most turbopumps are unsymmetricaily 
supported by manifolds, ducts, and brackets. This can 
be illustrated simply (Fig. 18) by representing a 
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Fig. 18. Simple Housing Model 


housing as a mass attached to ground with two radial 
degrees of freedom relative to the centerline as shown 
below. The equations of motion for Fig. 18 are: 



Assuming harmonic motion solution and appropriate 
derivations of: 



give the following solution: 
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Reformating the solution yields: 


(MX 4 c X 4- k ) 

XX XX 


(c X 4- k ) 
yx yx 


(c X 4 k ) 
xy xy 

(MX ^ 4 c X 4 k ) 

yy yy 


= o 


The stability can be investigated by the Routh- 
Hurwitz criterion (Ref. 4 ) as follows. The character- 
istic equation is : 


4 3 2 

a^X 4 a^X 4 a^X 4 a^X 4 a^ =0 


where assuming c^ - Cy X = 0, the coefficients for the 
above equation are: 


a o 

a i 

a 2 

a 3 

a 4 


= M(c + c ) 
xx yy 

= (k M + k M+c c ) 
yy xx xx yy 

= (c k + c k ) 

xx yy yy xx 

= (k k - k k ) 
xx yy xy yx 


The system is stable if the following quantities are 
greater than zero: 
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Substituting the above equations and expanding yields 

D = M(c + c ) 

1 xx yy 

D. = M^(c k + c k ) + M(c + c c 

2 xx xx yy yy XX yy xx yy 

= c c M^ (k - k + k k M^ (c + c 

3 xx yy xx yy xy yx xx yy 


+ Me c (c k + c c k + c c k 
xx yy xx yy xx yy yy xx yy xx 

+ c ,^k ) 

yy xx 


D 


4 


(k k 
xx yy 


k k ) 
xy yx 


3 


Note that, since k^, ky y , c xx , and Cyy are usually 
greater than zero, the possibility for a system insta- 
bility generally is controlled by the quantity. 
Therefore, in a symmetric case if k^ = -ky X , the sys- 
tem stability tends to decrease. However, if asym- 
metry exists then k xx 4 k^ and the first term on the 
righthand side of the equation is positive and the 
trend is for stability to increase. Since the rotor 
response is relative to the housing and the housing 
stability is increased, then the rotor /housing system 
must follow. This stability criterion cannot pre- 
dict the magnitude of the gain because in reality 
there are several coupling mechanisms that vary in 
some fashion with shaft speed. It was therefore nec- 
essary to construct a simple generic turbopump model 
where some of these ideas can be more easily evaluated. 


Generic Turbopump Model 

A model was constructed to simulate a typical 
turbopump with attention given to significant dynamic 
features . The finite element model consists of 10 
lumped mass and inertia locations for both the rotor 
and housing. Three large masses reside at appropriate 
locations in the rotor model to analytically represent 
two impellers and a turbine disk. A description of 
the model is provided in Fig. 19 . 
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1 2 flJl 4 5 [si 7 8 [ 9 I 10 

rotor f ot t mm' i <*) i w ♦ w i m 1 gnwr — ►* 


V BZ 


^BZ 


jll 12 13 14 15 16 17 18 19 f20 

housing hwymm m is * ♦ w* m vswju 

^ K sz ^ K sz 

SYMMETRIC IN Y AXIS 


MODEL: ROTOR SHAFT DIAMETER: OD = 3.0 INCHES, ID * 0.0 INCHES 

MATERIAL STEEL, E = 30 X 10 6 PSI 
SIMULATED LUMPED MASSES: POUNDS J13=10 

J16= 15 
J19= 20 


JOINT LENGTH - 3.0 INCHES 
ROTOR LENGTH * 27.0 INCHES 


HOUSING 

DIAMETER. OD * 15 INCHES, ID - 14 INCHES 
WEIGHT: 6/1 CASING/ROTOR WEIGHT RATIO 

k BZ = K BY = 500,000 LB/IN. 

K sz = 1 0 X 10 7 , 1/3 “ K sy /K sz 

Fig. 19. 20 Lumped Mass Turbopump Model 


Inputs to the analysis included mass and stiffness 
properties that were selected to be a generic repre- 
sentation of a turbopump. In addition, rotor/housing 
coupling coefficients for the HP0TP were arbitrarily 
selected as typical. These coefficients included the 
wear ring seal coupling on the first disk, the impeller- 
diffuser interaction on the second disk, and the 
turbine interstage seal coupling coefficients on the 
third disk. These coefficients are provided in Table 4. 

To further study rotor stability with a flexible 
housing, several dynamically important parameters 
were varied. The weight of the housing was varied 
from one to seven times the weight of the rotor. A 
housing/rotor weight ratio of approximately 6:1 is typ- 
ical of many existing turbopumps as shown in the data of 
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Table 4. Generic Turbopump Model Rotor/Housing 
Coupling Coefficients 

WEAR RING SEAL COUPLING ON FIRST DISK 


SPEED, 

nm 

*xx» 
LB/ IN. 

K xy» 
LB/ IN. 

’ c xx» 
LB-SEC /IN. 

1,000 

0 . 

0 . 

0. * ; 

10,000 

1.0 s 10 s 

1.5 x 10 4 

6.3 X 10 1 

20,000 

4.2 x 105 

5.6 x 10 4 

— 

25,000 

@.5 x 10 s 

8.6 x 10 4 

— 

30,000 

9,3. 1.10 s 

1.2 X 10 s 

2.1 x 10 2 

35,000 

1.2 x 10® 

1.6 x 10 5 

— 

40,000 

1.6 x 10® 

2.0 x 10 s 


45,000 

2.1 x 10® 

2.7 X 10 s 

3.2 x lO 2 


IMPELLER - DIFFUSER INTERACTION 
ON SECOND DISK 


SPEED, 

^xx* 
LB /IN, 

Kxy, 
LB /IN. 

c xx* 1 

LB-SEC /IN. 

1,000 

-8.9 x Id 1 

-4,0 x Id 1 

0.1 

10,000 

-0.® x I© 3 

-4.0 x 10 3 

1.0 | 

20,000 

-3.5 x I© 4 

-1.6 x I® 4 

2.0 i 

2$, TOO 

-5.5 x 10 4 

-2.5 x 10 4 

2.5 

30,000 

-8.0 x 10 4 

-3.8 x I® 4 

3.0 

35,000 

-1.1 x 10® 

-4.9 x 1§ 4 

3.5 

40,000 

-1.4 x 10 s 

-8.4 x Id 4 

4,0 

45,000 

-1.8 x I© 5 

-8.1 x I© 4 

45 j 


TURBINE INTERSTAGE SEAL COUPLING 
ON THIRD DISK 


j SPEED, 

nm 

LB /IN. 

lb /In. 

LB-SEC / IN. 

1,000 

0 . 

0 . 

0 . 

10,000 

8.7 x 10 3 

1,9 x 10 3 

2.3 

20,000 

3.5 x 10 4 

7.9 x 10 3 

6.1 ! 

25,000 

4.7 x 10 4 

1.3 x I® 4 

8.8 | 

30,000 

6.8 x I® 4 

2.1 x IQ 4 

12.8 ; 

35,000 

9.5 x 10 4 

3.1 x 10 4 

17.3 1 

40,000 

1.3 x 10 5 

4.6 x 10 4 

23.0 

45,000 

I-,.. J 

1.6 a I® 5 

6.4 x I® 4 

29.6 
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Table 5. Since the significant dynamic difference be- 
tween a rigid and flexible housing is asymmetric modes 
and with insight provided by Ref. 5 , the effects of 
support asymmetry was determined by varying the support 


Table 5. Housing/Rotor Weight Ratio 
of Existing Turbopumps 


MODEL 

ROTOR WEIGHT, 
POUNDS 

HOUSING WEIGHT, 
POUNDS 

HOUSING /ROTOR 
WEIGHT RATIO 

MK38-0 (104%) 

81. 

480. 

6/1 

MK38-0 (109%) 

81. 

480. 

6/1 

MK38-F 

128.5 

649. 

5/1 

MK48-F 

8.25 

87.75 


MK49-F 

6.1 

48. 

8/1 


stiffness ratio in orthogonal directions. Figure 20 
shows the stability threshold speed plotted as a func- 
tion of the support asymmetry ratio for a family of 
housing/rotor weight ratios. Results show that there 
is an optimum value of asymmetry that produces the 



Fig. 20. 20 Lumped Mass Model (Stability 

Threshold vs Support Asymmetry 
Ratio) 
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maximum stability threshold speed for a given housing/ 
rotor weight ratio. This optimum asymmetry ratio varies 
as the housing/rotor weight ratio is increased. The 
larger the housing/rotor weight ratio, the smaller the 
optimum asymmetry ratio to produce a maximum threshold 
speed. For example, an asymmetry ratio of 3:1 produces 
the highest threshold speed for a large housing/rotor 
weight ratio, whereas an asymmetry ratio of 10:1 re- 
quires a low weight ratio to produce the highest thres- 
hold speed. 

Finally, housing modal damping was varied to de- 
termine its effect on rotor stability. The modal damp- 
ing (percentage of critical damping) supplied by the 
housing structure was varied with increasing asymmetric 
support ratios. Figure 21 shows that as the housing 
modal damping is increased, the stability threshold 
speed remains nearly constant. This result does not 
change with the housing/rotor weight ratio. 

38,000 


37,500 


oc 37,000 

d 


ui 

| 36.500 
V- 
> 
t 

-j 

< 36.000 

H 

to 

35,500 


35,000 

0 .010 .020 .030 .040 .050 .080 .070 .080 .090 .100 .110 

MODAL DAMPING 

Fig. 21. Hodal Damping vs Stability Threshold 
(20 Lumped Mass Model) 
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In an effort to compare the previously observed 
trends of the six turbopumps to this generic model , 
stability analyses using both a rigid and flexible 
housing assumption were performed. Figure 22 shows a 
stability map for the simple model. The flexible 
housing approaches an unstable condition at a higher 
speed than the rotor in a rigid housing. This result 
emulates the turbopump survey performed earlier. The 
analysis was performed with the typical 6:1 housing/ 
rotor weight ratio and a support asymmetry ratio of 
3:1. 



Fig. 22. Generic Turbopump Model 
Stability Map 


Conclusions and Recommended Design Practice 

A rotor stability analysis with a rigid housing 
assumption predicts an instability at a lower speed 
compared to a similar flexible housing case. Conse- 
quently, a rigid housing assumption is more conserva- 
tive for design than the flexible model. Simple 
analysis approaches indicate that a flexible housing 
acts as an inefficient energy absorber and dissi- 
pates only a small portion of the destabilizing 
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forces in the rotor /housing system. The significant 
stability improvement associated with a housing is 
attributed to asymmetric supports. In a simple generic 
turbopump model, it was determined that the asymmetric 
support stiffness ratio grounding the housing should be 
greater than a value of 2:1. This simple model also 
indicated the optimum housing asymmetric support ratio 
is dependent on the housing/rotor weight ratio. Review 
of the weight data for the stable turbopumps included 
in the survey and the weight parametric study indicate 
most well behaved pumps have or should have at least a 
5:1 housing to rotor weight ratio. Housing structural 
damping variations indicate no significant effect on 
the rotor stability. 

For future designs , these results indicate some 
general rules to act as guidelines . These include: 
asymmetric housing supports of greater than a 2:1 
ratio and heavier than a 5:1 housing to rotor ratio. 
Also, preliminary analyses should predict conservative 
stability results assuming a rigid housing depending 
on the accuracy of the coefficients that couple the 
rotor to the housing. Based on analytic results, in 
many cases complex housing models are not necessary 
for stability analyses. As a result, time and money 
are saved in the development of rotordynamic models . 

Housing models would be required if detailed re- 
sponse analyses are necessary, however , the variation 
in response with unknown rotor unbalance distributions 
many times makes preliminary response work unreliable. 
In practice, a high-speed rotor balance in the housing 
is recommended to obtain minimal response. 

Nomenclature 

A,B = constants, determined by initial conditions 
c = damping, lb-sec/in. 

F - excitation force, pounds 

k « stiffness, lb /in. 

2 

m,M = mass, lb-sec /in. 
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arbitrary constants 
ro tor displacement , inches 
rotor velocity, in. /sec 
housing displacement, inches 
housing velocity, in. /sec 
time, seconds 

work per cycle, lb- in. /cycle 
spin speed, rad/sec 
displacement, inches 
velocity, in. / sec 
acceleration, in. / sec 
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